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^ . Doubly charged scalars, predicted in many models having exotic Higgs represen- 
I tations, can in general have lepton-nuniber violating (LFV) couplings. The basis 



of most searches for this charged scalar has been to look for its direct production 
and its subsequent decay to like-sign final state leptons. In this work we show 
that by using an associated monoenergetic final state photon seen at a future lin- 
ear e~e~ collider, we can have a clear and distinct signature for a doubly-charged 
resonance and also determine its mass rather precisely We also estimate the 
strength of the AL = 2 coupling which can be probed in this way at y/s = 1 
TeV, as a function of the recoil mass of the doubly-charged scalar. 



The absence of the Higgs boson from the cupboard containing particle trophies discovered 
by experimentahsts, still leaves a scope of speculation as to what would eventually be the 
possible structure for SU{2) x t/(l)-breaking. Thus scenarios with extended Higgs sectors, 
ranging from ones with two or more doublets to those with other representations of SU(2), 
are often considered. Doubly charged scalars arise in a number of such scenarios [1,2]. The 
most common models to accommodate such scalars are those with triplet Higgs. Triplets 
can be made part of the electroweak symmetry breaking sector in purely phenomenological 
studies [3-5] . On the other hand, they may be indispensable in some special theories such as 
the simplest versions of Little Higgs models [6], where the Higgs is envisioned as a pseudo- 
goldstone boson, and the triplets are required to cancel the quadratic divergence to the Higgs 
mass. 

An added feature often associated with doubly- charged Higgs is the possibility of lepton- 
number violation. This basically consists in AL = 2 couplings with leptons of the form 



Cy = ihij%LCT2^-^ jL + h.c. (1) 

where i,j = e,/i, r are generation indices, the \E''s are the two-component left-handed lepton 
fields, and $ is the triplet with Y — 2 weak hypercharge and is given by the 2x2 matrix of 
the scalar fields: 

^ (j>+/^/2 0++ ^ 

This leads to mass terms for neutrinos [1,7] once the neutral component 0° acquires a 
vacuum expectation value (vev): 

~ K,v' (2) 

v' being the triplet vev. Since constraints on the p-parameter [8] puts strong limits on the 
the triplet vev [9] in general, this immediately translates to limits on the L- violation Yukawa 
couplings from the expected ranges of neutrino masses [10]. Such limits usually constrain the 
collider signals for doubly-charged scalars sought through AL = 2 interactions. Of course, 
there are models where the limits from the p-parameter can be avoided [11] by postulating 
real as well as complex triplets at the same time, and assuming a custodial symmetry relating 
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their vev [3] . Although such a custodial symmetry has been found to be preserved in higher- 
order corrections from the scalar potential [12], its stability against corrections via gauge 
coupling is not obvious. Thus it is safe to abide by the constraints on both the triplet vev 
(0°) and the quantities hij{4>^) in an analysis related to collider signals. 

In this paper we point out the usefulness of looking for doubly-charged scalars in an e'e" 
collider, in the radiative production channel. A linear collider with, say, y/s — 1 TeV can 
operate for part of its run-time in the electron mode, where certain signals can be remarkably 
free from backgrounds. Also, this mode is perhaps ideal for exploring scenarios with AL = 2 
couplings. As for doubly-charged scalars, their resonant production in e~e~ as well as 
have been already studied [11, 13]. However, resonant production of the (f) requires one to 
know its mass with reasonable accuracy to start with, and tune the center-of-mass energy 
of the colliding electrons accordingly. In addition, precise identification of a doubly- charged 
resonance will also depend on its decay products [14], which depend on the parameters of 
the L-violating sector. In general, one can have the decays 

• 0— — > w-w- 

• 0— — > ri- 

• 0~ — > w-(t)- 

• (f) — > 0^0^ 

Of these, the third mode, if kinematically possibly, is dominant as it is driven by gauge 
coupling. However, a degeneracy among the triplet components is often a consequence of 
theories, albeit in a model-dependent fashion. If we thus neglect the last two channels hsted 
above, we still have the W^W~ and /~/~ channels, of which the first is controlled by the 
triplet vev v' and the second, by the coupling hu. When the first mode is dominant, it 
requires careful analysis of the W-decay products in order to isolate signatures of resonant 
production. Furthermore, the analysis becomes complicated in case the can decay into 
a 0~ or a pair of them. 

It is thus desirable to have supplementary channels in mind while looking for doubly- 
charged scalars. With this in view, we have calculated the rates for the process 



concentrating on the hard single photon in the final state. This photon will be monochro- 
matic if a doubly-charged resonance is produced, irrespective of what it decays into. Further- 
more, one is no more required to tune the electron-electron center-of-mass energy at a fixed 
value [15]. And finally, one can use such a study of monochromatic single photons in e~e~ 
collision to include the search for doubly-charged objects other than the scalar discussed 
above, an example being bilepton resonances [16]. 

For our analysis, taking the radiative production of the scalar (f) as the benchmark 
process, we concentrate only on the flavor diagonal coupling hee- Others, especially the non- 
diagonal couplings, are subject to very stringent bounds from rare decay processes of /x^ and 

leptons. The rare decay studies [17] provide bounds on the product of couplings as: 

he^,Ke < 3.2 X 10"^^ GeV-^Mjii (3) 
/ieA^ < 2.0 X 10-^° GeV-2M|±± (4) 

The relevant bounds in our case are the following upper bounds, which come from Bhabha 
scattering [18]: 

hl^ ~ 6.0 X 10-^ GeV-^M|±± (5) 

and from {g — 2)^ measurements [19]: 

~ 2.5 X 10-^ GeV-^M|±± (6) 

Much stringent upper bounds exist from the measurements of muonium-antimuonium tran- 
sition [20] in the form of the product of the couplings: 

KeKi, ~ 2.0 X 10-^ GeV-'M|±± (7) 

The above bounds allow for small value of doubly-charged scalar mass with small coupling 

constant. In our numerical estimate, we have chosen the coupling strength to be hee = 0.1 

which is consistent with the hmits obtained from Bhabha scattering [21,22] and also respects 

the most stringent bounds coming from muonium-antimuonium conversion results which for 

flavor diagonal coupling is h < 0.44 TeV~^ at 90% C.L. The latter bound however, can 

be relaxed in different new physics models [23]. 

As has been mentioned already, on-shell radiative production of a doubly-charged scalar 

gives an almost monochromatic photon of energy 

s - ML- 



which stands out against the continuum background of the standard model (SM). In the 
discussions to follow based on our work, we show how, by simply tagging on the isolated 
photon in the final state without bothering about the the other associated particles, and 
looking at the line spectrum superposed on the continuum background, leads to clear signals 
for the production of a doubly-charged scalar. 

Before we present the results of our analysis, it may be noted that there exist studies on 
the doubly-charged scalar in the context of hadron colliders [24] and in different modes of 
operation of linear colliders other than e~e~ collision [13,25] and e+e^ mode [26], such as the 
67 and 77 modes [27]. Although these works point out interesting ways of producing doubly- 
charged scalars either singly or in pair, the latter are somewhat restrictive in mass reach, 
considering the fact that a back-scattered photon carries a fraction of the initial electron 
beam energy. Moreover, there is an unavoidable kinematic suppression in pair-production in 
the 77 channel, in spite of the spectacular enhancement of rates due to the electric charge(s). 
Also, another possible channel for associated production of a doubly charged Higgs in e~e~ 
collisions has been looked at in ref [28]. In this work we show that we can probe the full 
energy reach of the collider in its fundamental mode of running, and is only restricted by the 
phase-space restrictions due to the kinematic cuts used for the selection of events. Hardness 
and transversality cuts on the photon in our final state are presumed to avoid any confusion 
with initial state radiation (ISR) or beamstrahlung photons. 

Y 

Figure 1: Feynman graphs corresponding to single production of a doubly-charged scalar along 
with an associated photon in e~e~ collisions. 

The Feynman diagrams contributing to the associated photon process are shown in Fig- 
ure m It is worth mentioning here that since they involve clashing fermion lines, one has to 
be careful while writing down the Feynman amplitudes and carefully handle the charge con- 
jugation operator appearing in the Feynman rules for the eeH vertices. Using the properties 
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of the charge conjugation operator C, it is a matter of simple algebra to write down the 
matrix amplitude squares for the different graphs in Figure H The above process is also a 
clear indicator of a AL = 2 process. We look at this process in the context of a ^/s = 1 TeV 
linear e~e~ collider. As discussed earlier, we consider decay of the doubly-charged scalar to 
like signed leptons only and that too of the same flavor. We assume that the ha couplings 
are of equal strength and so the branching ratio BR(0 ^7^7) = 1/3 for i = 1,2,3. In 
this work we have consistently chosen the coupling strength to be hee = 0.1 and we calculate 
the decay width of the doubly-charged scalar assuming the decays to leptons only. In other 
words, we assume the triplet vev to be very small, so it hardly contributes to the decay 
mode of (p W~W~ which is directly proportional to the triplet vev squared (f'^) and 
we have also neglected the other possible decay modes of cf) as pointed out earlier. We find 
that the total decay width obtained is very miniscule (~ 1.2 GeV) for a 1 TeV scalar mass 
when compared to the machine energy. So we can use the narrow-width approximation and 
consider on-shell production of the doubly-charged scalar and its subsequent decay to e~e~. 

The major SM background that contributes to the above process is the radiative Moller 
scattering process: 

e~ -|- e~ — >• 7 -|- e~ -|- e~ 

which, although a continuum background, could prima facie be large enough to wash away 
the monochromatic peak. The event selection criteria, therefore, are largely aimed at sup- 
pressing this continuum background. 

We impose the following set of cuts. Since we do not want the final state particles to be 
too close to the beam pipe, one needs to have a rapidity cut on the final state particles: 

|r/(e")|<3.0 and |r/(7)| < 2.5 

We also demand a hard photon in the final state, which is the main focus of our work. This 
is obtained by imposing a cut on the minimum photon energy: 

E(7) > 20 GeV 

A minimum energy is also demanded for both the final state electrons: 

E{e~) > 5 GeV 
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Figure 2: Differential cross-sections for (a) signal with \rje\ < 3.0, (b) signal+background and SM 
background with \r]e\ < 3.0 and (c) signal+background and SM background with |r/e| < 1-5, against 
electron rapidity rje ■ 



The above criteria help us in suppressing the continuum background to a considerable extent. 
Finally we would like the detectors to register and resolve events for the different particles 
and hence all the final state particles should be well separated in space and satisfy: 

5R > 0.2 

where [SR)"^ = (A0)^ + {Ar])"^ with Arj and A0 respectively denoting the separation in 
rapidity and azimuthal angle for the pair of particles under consideration. Using the above 
cuts we make an estimate of the SM background and the signal. The background has been 
generated using the MadEvent Monte Carlo generator [29] . It is found that the background 
is quite large compared to the signal with a cross-section of about 3216 fb. In Figure El we 
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show the distribution for the differential cross-section as a function of the electron rapidity. 
It is worth noticing that the SM background is symmetric in the binwise distribution for 
?7(e~) of the final state electron and is peaked away from the central region of the rapidity 
distribution (Fig|2fb)), which means that they are more in the forward direction. This is 
expected due to the strong t-channel radiative contribution to Moller scattering. On the same 
curve we also show the distribution for the signal+background for two different values of the 
doubly- charged scalar mass, — = 300 GeV and 750 GeV. In contrast to the background, 
we find that the signal is peaked at the central region, which is clear from Fig|21^a) where 
we have plotted the signal alone and this peaking becomes more pronounced as the mass 

3 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I i| 
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Figure 3: Illustrating the ratio of signal+bacliground with the SM bacliground. The broken line 
(blue) corresponds to the cut |r/e| < 3.0 and the solid line (red) corresponds to the cut \r]e\ < 1.5 

of the increases. Infact for the low mass there does not seem to be much difference 
in the rapidity distribution but as the mass is increased to a higher value, the difference 
shows up with the signal peaked centrally. This can be easily understood if we look at the 
kinematics of production. The higher the mass of the (p , the less boost it will have and 
hence the decay products will come out back-to-back. Using this as a cue, we can actually 
implement a more stronger cut on the rapidity of the electron which will throw the signal 
into prominence. We find that a cut of 

|r/(e-)|<1.5 
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is good enough to kill the background by about 80% while the signal goes down only by 
about 45% for low masses (~ 300 GeV) and by about 10% for large masses (~ 900 GeV) for 
the doubly-charged scalar. With the kind of luminosity expected at future linear colliders 
we expect that the event rates would be considerable even if low mass state is realized and 
although the signal might drop by 40% — 50%, it will still stand out against the much reduced 
background due to this cut. To highlight this we implement the cut and plot the distribution 
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Table 1: Cross-sections for signal and the SM background corresponding to the two choices of 
rapidity cut on the final state electron. 

in Fig|21^c). It is clearly visible how the signal is enhanced compared to the background with 
this cut. We have demonstrated the signal+background for three different masses for the 
doubly- charged scalar, — = 300, 600, 900 GeV. We also show the significance of the cut 
in Fig El by plotting the ratio of the cross-sections of the signal+background (as+B) with the 
SM background (cb), for both the choices of the cut on electron rapidity, as a function of the 
doubly-charged scalar mass. The graph clearly highlights the enhancement in the signal to 
background ratio and reflects the increase in the cross-section as we go higher in scalar mass. 
In Table 1, we list the cross-sections (rounded off to the nearest integer) for different masses 
corresponding to the respective cuts imposed on the electron rapidity for a more quantitative 
outlook. 

Next, we focus on the main trigger, viz. the photon. In Fig H^a) we show the distri- 
bution of the photon energy, where we have superposed the differential cross-section for 
signal+background in each bin over the SM background. A pronounced peak can be seen in 
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Figure 4: Differential cross-sections against (a) pfioton energy and (h) invariant mass of electron 
pair Mee- The dash-dot-dash (green) line corresponds to — = 300 GeV, dotted (blue) line 
corresponds to — = 600 GeV and the dashed (red) line corresponds to — = 900 GeV 
respectively. The binsize is chosen to be 5 GeV in (a) and 10 GeV in (b). 



the photon energy distribution, due to the monochromaticity of the photon, corresponding 
to the recoil energy against the scalar resonance through the relation of Eq|S] To make 
our analysis realistic, we have smeared the photon energy by a Gaussian function whose 
half-width is guided by the resolution of the electromagnetic calorimeter [30,31]: 

AE 14% 



E 

Moreover, we have fully incorporated the effects of ISR [32] which often results in substantial 
broadening of the peak, due to the spread in the effective center of mass energy available for 
our process. We have used CompHEP [33] to include ISR effects for the SM background. We 
show the resulting peak for three choices of scalar mass (300, 600, 900 GeV). Alternatively, 
in FiglU^b), we also show the invariant mass distribution of the ee pair for the above choice 
of parameters and as expected the distribution peaks corresponding to the mass of scalar. 
Here the half-width of the Gaussian function used for smear is [30,31]: 

AE 15% 

In Fig El we plot the energy distribution of the photon once again. But here we assume 
that we do not have prior knowledge of the decay products of the doubly-charged scalar. In 
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Figure 5: Differential cross-sections against photon energy when cf) X (anything). The 

dash-dot-dash (green) hne corresponds to M^~~ = 300 GeV, dotted (blue) line corresponds to 
— = 600 GeV and the dashed (red) line corresponds to — = 900 GeV respectively. The 
binsize is 5 GeV. 

other words, we only look at the final state hard transverse photon in 

e~e~ — i> 7 + — > 7 + X 

The distribution again shows peaks corresponding to the recoil against the massive scalars, 
irrespective of the knowledge of the decay products of the scalar. In fact our signal here 
receives a relative boost as it is not suppressed by considering any further decay since the 
BR{(f) X) = 100%. Through FigEJ and Fig El we have shown that a single associated 
photon will show peaks in its energy distribution over the continuum background of SM, if a 
doubly-charged scalar is produced with mass — < y/s. We only demand that no electron 
in the final state can have a rapidity whose absolute value is greater than 1.5. The fact 
that looking at a single photon against the backdrop of a continuum background makes it 
possible to identify a LFV (AL = 2) process in a model independent way, makes this signal 
worth studying at a future e~e~ collider and running the linear collider in this mode. 

Since the rates for the signal depend directly on the eeH coupling squared, we can make 
an estimate of the strength of the coupling which can give substantial rates for identification 
of the peaks in the photon energy distribution. We do this for the case when X 
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Figure 6: Illustrating the reach of the coupling constant at which the resonances in the distri- 
bution can be identified over the fluctuations in the SM background. The assumed luminosity is 



and in Fig IHl We show the strength of the couphng for which the peaks would stand out 
against the fluctuations in the SM background at 2a, 4a and 6a level. In our analysis we 
have assumed a luminosity of C = 100 fb^^, easily achievable at future linear colliders. We 
have also done this analysis based on the smear in photon energies due to finite resolution 
of the detector as well as ISR effects, as mentioned earlier. A bin of width 10 GeV about 
the signal peak in photon energy has been identified for each — . We then look at the 
fluctuations in the SM background in that bin and compare the corresponding rate for the 
signal. This procedure has been repeated for different coupling strengths. The fact that 
we are not looking at any specific final state arising from decay improves the reach of 
this search channel. Our analysis suggests that the peaks in the photon distributions will be 
distinguishable for coupling strengths as low as 0.006(0.010) for scalar mass of 300 GeV at 
2a{6a) level and to about 0.0034(0.006) for scalar mass of 900 GeV at 2a{6a). This estimate 
far overwhelms the simple method of comparing total rates of signal and background and 
restricting the coupling strength in the parameter space. It is also worth mentioning that 
if one is to assign a certain order of detection efficiency e, with the final states then the 
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above reach atmost scales by e~^/^. However, if a direct resonance is excited then that would 
invariably translate into a much stronger probe of the coupling strength [11]. Nonetheless, 
our analysis is not dependent on the tuning of the ^/s of the machine to hit a resonance and 
hence serves as a more robust proposition. For luminosity higher than what we have used, 
this reach can be further enhanced. 

To summarise, the cleanliness of central photon detection at a high energy linear collider 
can be very helpful in identifying a doubly- charged scalar. While it is true that final states 
such as {fi~fi~) can be completely background-free, the peaks in the hard photon energy can 
be helpful in two ways. First, one does not need to tune the two electron beams, and can 
therefore work without a prior knowledge of the (p mass. Secondly, this method is shown 
to work even if the dominantly decays into states that are not clean enough for the 
resonance to be identified. Thus, as soon as one succeeds in reducing the SM backgrounds 
(using, for example, the electron rapidity cut), one can clearly see AL = 2 interactions, just 
by looking at the accompanying hard photon. Not only doubly-charged scalars but also more 
exotic resonances such as bileptons are amenable to detection in this manner. 
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